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Abstract 

The feasibility of an adaptive control system designed to 
alleviate helicopter gust-induced vibration has been analytically 
investigated for an articulated rotor system. This control system is 
based on discrete optimal control theory, and is composed of a set of 
measurements (oscillatory hub forces and moments), an identification 
system using a Kalman filter, a control system based on the minimiza- 
tion of the quadratic performance function, and a simulation system 
of the helicopter rotor. The gust models are step and sinusoidal 
vertical gusts. Control inputs are selected at the gust frequency 
(wq) , subharmonic frequency (wq - fi) , and superharmonic frequency 
(aiQ + , and are superimposed on the basic collective and cyclic 

control inputs. The response to be reduced is selected to be that at 
the gust frequency because this is the dominant response compared 
with sub- and superharmonics. Numerical calculations show that the 
adaptive blade pitch control algorithm satisfactorily alleviates the 
hub gust response. Almost 100% reduction of the perturbation thrust 
response to a step gust and more than 50% reduction to a sinusoidal 
gust are achieved in the numerical simulations. 


Nomenclature 

a 

b 

Cr » Gljl , Cy 

•^MX’^MY’^MZ 

c 

E( ) 

Fx.Fy,Fz 

8 


lift slope 
semichord, c/2 

horizontal, thrust, and side force coefficients 

rolling, pitching, and yawing moment coefficients 

chord length 

expectation of ( ) 

hub forces 

gust frequency, Hz 

gust velocity vector 


*NRC Associate. 
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^3 

moment of inertia of a blade about flapping hinge 
= J (r - rg)^m dr 


^X’My*Mz 

hub moments 


“3 

m 

mass moment of blade = J r(r - rg)m dr 

’^3 

blade section mass 


Q 

variance of the process noise 


r 

variance of the measurement noise 


’'3 

position of flapping hinge 


t 

time 


Ug.Vq.Wg 

horizontal, lateral, and vertical components of gust 
velocity 


V 

measurement noise 


w 

gust amplitude, m/sec 


z 

system response measurements 


A 

perturbation from a steady value, or small Increment 


P 

advance ratio 


Q 

rotor rotational speed, rad/ sec 


(Hq 

gust angular velocity, rad/sec 


Subscripts: 



C 

control, or cosine element 


G 

gust 

' 

n 

n time-cycle 


S 

sine element 


0 

initial value, or amplitude, or uncontrolled value 


Superscripts: 



(“) 

trimmed value 


(•) 

d( )/dt 
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estimate of ( ) 


T 

( ) transpose of ( ) 


1. Introduction 


Helicopters fly close to the ground where the motion of the 
atmosphere may be thought of as turbulent flow. To avoid structural 
vibrations and unfavorable dynamic characteristics in flying and 
riding qualities caused by gusty winds, it is necessary to analyze 
the gust response of a rotary wing aircraft and develop a scheme to 
alleviate such responses. 

Studies of gust responses of a helicopter rotor were reported 
in Refs. [l]-[5]; only Refs. [3] and [5] dealt with rotor gust 
response experiments. Experimental data for vertical gusts show 
that the thrust response was more sensitive than other responses of 
the rotor. For the sinusoidal gust with frequency (uiq) , the thrust 
response characteristics have not only a (wq) component, but also 
(o)Q + nQ) (n = ±1, ±2, . . .) components in the fluctuation. In this 
thrust response, the (ioq) element of the oscillatory characteristics 
is dominant in the comparatively low gust frequency range (cdq « Q) . 

Attempts to alleviate the gust-induced vibration have been 
described in Refs. [6]-[10]. Among them, Briczinski and Cooper [6] 
have investigated the effect of a rotor/vehicle state feedback system 
on the handling qualities of a helicopter, specifically characteris- 
tics concerned with gust response. They found that the feedback 
scheme of the rotor tip-path-plane or body state was very useful as 
a means of gust suppression. Johnson [7] studied the performance of 
an optimal control system applied to proprotor/wing gust response. 
Significant and simultaneous reduction in the rotor and wing responses 
was achieved. Zwicke et al. [8] investigated the performance of an 
optimal sampled-data feedback system on the gust response. They also 
studied a suboptimal feedback system derived from the above control 
system; a significant reduction in the gust response using this sub- 
optimal feedback system was achieved. Ham and Mcklllip [9] developed 
an individual blade control (IBC) scheme for gust response allevia- 
tion. Saito et al. [10] also studied a simple feedback system to 
alleviate rotor gust response; in their control scheme. Individual 
blade pitch angle control, based on scheduled feedback gains derived 
from analytical calculations, is used. 

There has been great progress in active control vibration 
reduction techniques for helicopters (theoretically and experimen- 
tally) in the past decade (Refs. [11]-[15]). In these vibration 
reduction systems, control schemes known as multicyclic control 
(Ref. [11]) or higher harmonic control (Ref. [13]) have been applied 
to the reduction of inherent vibratory responses of a helicopter. 

Pitch Inputs at harmonics of the rotor rotational speed are used. 
Typically the helicopter is represented by a linear, quasi-static 
frequency domain model. The relationship between control inputs and 
outputs (which can include vibrations, loads, and rotor performance 
parameters) are modeled by a transfer function matrix. Theoretical 
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and experimental results show that the vibratory level of a helicopter 
can be significantly reduced by using controllers in which the trans- 
fer function matrix is updated by a Kalman filtering scheme. These 
vibration reduction systems are reviewed in detail in Ref. [15]. 

In this paper, the feasibility of an adaptive control system 
designed to alleviate gust-induced vibration has been analytically 
investigated. This control system is based on a discrete optimal 
control theory, and is composed of a set of measurements (oscil- 
latory hub forces and moments), an identification system using a 
Kalman filter, a control system based on the quadratic performance 
function, and a simulation system of the helicopter rotor. The 
gust models are step and sinusoidal vertical gusts. The local 
momentum theory (LMT) (Ref. [16]) is used to calculate time-wise 
vibratory airloads and moments at the hub position of the rotor. An 
H-34 rotor with four articulated blades is used. The blades have 
full flap, lead-lag and torsion elastic deflection. Fuselage motion 
is not considered. 


2. Adaptive Blade Pitch Control 

An adaptive blade pitch control designed to alleviate gust- 
induced vibration is analytically investigated. For the gust response 
of a helicopter, additional vibratory responses, such as gust harmonic 
(ojq) , subharmonic (ioq - fi) , and superhamonic (o)q + f^) , etc., appear 
in the response. Gust harmonic response is dominant in the compara- 
tively low gust frequency range (Refs. [3], [5]). Hence, this 
response is selected as the response to be reduced. In this study, 
the helicopter’s gust response is characterized by time-dependent hub 
forces and moments. 

Helicopter Model 

In this study, the helicopter is represented by a linear, 
quasi-static frequency-domain model relating the output {z} to the 
input {0}. Here {z} is an amplitude vector of the gust-induced 
vibration harmonics in the nonrotating frame. The input {0} is at 
the gust frequency. It is assumed that the gust is sinusoidal, and 
that the gust frequency is known a priori. When a rotor penetrates 
into the gust, the hub oscillatory forces and moments can be consid- 
ered to be composed of gust and control components: 

{z} = {zq} + {z(;;} (1) 

Equation (1) can be rewritten by using the expression of the rotor 
Impedance [Tq] and rotor transfer function [T^] as follows: 

{zg} = [TG]{g> 

{zcJ = [Tc]{0} 

That is, Eq. (1) becomes 
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{z} = [TG]{g> + [TcHe} 


( 2 ) 


where 

{z} = (m X 1) vector 

[Tq] = (m X 3) matrix 

{g} = (3 X 1) vector 

[T^] = (m X n) matrix 

{0} = (n X 1) vector 

Once the hub response characteristics to a gust are determined, [Tq] 
becomes a constant matrix. Therefore, we may denote the first term 
of Eq. (2) as {zq}. Then Eq. (2) is 

{z} = {zj + [Tc]{0} (4) 

This form resembles the global model formulation of helicopter vibra- 
tion (Ref. [15]). The following "local model" can be also taken into 
account: 

{Az} = [TcUAe} (5) 

where z means z^ - z^-i and A6 means 0^ - ©n-i • 

Oscillatory hub forces and moments at the gust frequency in the 
nonrotating frame are chosen as the measurements {z}: 

{z} = (F 2 c, F 2 x, Fxc. Fxs. Myc* Mys, Fyc, Fyg, Mxc. Mxs. MgC. Mgs) 

( 6 ) 

where the subscripts C and S denote the cosine and sine components 
of the hub reaction at the gust frequency wq. Control inputs are 
selected at the gust frequency (o)q) , subharmonic (wq - 0) , and super- 
harmonic frequency (ojq + in the nonrotating frame as follows: 

A0 = 0^ cos(wQt) + 02 sin(wQt) + 0g cos(wq - 0)t 

+ 0^ sin(u)Q - fi)t + 0g cos(o)q + fi)t + 0g sin(o)Q + 0)t (7) 

Hence the gust-control vector {0} has the components 

{0} = (01, 02, 03, ©4, ©5, eg) (8) 

These control inputs are then superimposed on the blade trim pitch 
inputs. 
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Identification 

Three cases are distinguished for the helicopter model, depend- 
ing on the identification approach: 

i. Identify {zq} only 

ii. Identify [Tq] only 

iii. Identify {Zq} and [T(j] 

For the gust response, the uncontrolled response {zq} is generally 
time-variant; the matrix [Tg] depends on the operating flight condi- 
tions. Hence it is necessary that a transfer function [Tg] be iden- 
tified simultaneously with the uncontrolled response {Zg}. Case ill 
is considered in this investigation. For the local model, it is 
necessary that the [T^] matrix be identified for each time-cycle. 

The Kalman filtering technique (Ref. [17]) is applied for identifica- 
tion. For the global model, Eq. (4) can be rearranged as follows: 

Zn = Zq + Tq6j^ = [Tq, z] (6jj» 1) 

= T • 0n (9) 

where 

T = m X (n + 1) matrix 
9n = (n + 1 ) X 1 vector 

The subscript n denotes the time-step at t = n At. In this study, 
it is assumed that there is no noise in the measurement of 9^^. The 
identification algorithm can be derived by considering the jth mea- 
surement as 


jn 


= 9^t. 
n jn 


+ V 


jn 


( 10 ) 


where tj^ is the jth row of T and vjn is measurement noise, 
which has zero mean E(Vn) = 0, variance ECv^Vm) = I’n'^nm* ^ 
Gaussian probability distribution. Here is the Kronecker delta 

function and the subscript j will be omitted to simplify the nota- 
tion. The variation of the parameters is modeled as a random process: 


n+i 


= t„ + u 


n 


( 11 ) 


where u^j is a random variable with zero mean E(Uj^) = 0, variance 
E(uj^Ujjj) = Qn^nm* ^ Gaussian probability distribution. The minimum 
error -variance estimate of t^j is then obtained from the Kalman 
filter (Ref. [17]). 


*-n ^n-i + ^i^^n ~ ®n^n-i) 


( 12 ) 
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where 


(13a) 

(13b) 

(13c) 


+ Q„.. 

Pn " «n - “n*n*iV<'n + *i«n*n) 
'‘n = Vn/<Pn + ‘SVn) 


Here is the variance of the error in the estimate of t^ before 

the measurement, and Pn is the variance after the measurement. The 
parameter is the Kalman gain vector. To simplify the calcula- 

tion, it will be assumed that and rjj have the same time variation 
for all measurements, and that Qn» proportional to 

the same function fj: rjj^ = Ijr^; Qjn = ^jQnJ ^jo “ Then it 

follows that = f ^jn “ ^j^n’ that the gain kj^ is 

the same for all the measurements. With the same gains, the rows can 
be combined to form 


IT., = ITr, 

n n-i ' n 


n-i n' n 


(14) 


For the local model, Eq. (14) takes the following form: 

^n ~ ^n-1 f^^n ~ ^n-i^ “ ^n-i ^®n ~ ®n-i^^^n 

where 

^n “ ^n-i Qn-i (16a) 

Pn ■ - «n *«n (16b) 

■'ll ■ “n 66n/(r„ + ae^ a«n) 

= 9n ' «n-i d6'« 

Controller 


The control algorithm is based on the minimization of a per- 
formance index J that is a quadratic function of the input and 
output variables. The quadratic performance function used here is 

J = ^X^n + ^X^n + ^0X0 A9n (17) 

where Wg, W 0 , and W ^0 are weighting matrices, which are assumed to 
be diagonal, and all harmonics of a particular quantity have the same 
value. Then J is a weighted sum of the mean squares of the gust 
response and control. The matrix W 0 constrains the amplitude of the 
control and W ^0 constrains the rate of change of the control. 

For the deterministic controller, the control required to 
alleviate the helicopter vibration is given by substituting for 
in the performance function J using the helicopter model and then 
solving for 9^^ that minimizes J. For the global model (Eq. (4)) 
the solution can be obtained as follows: 
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( 18 ) 


where 


C . -dtJw^ 

^A0 “ 

0 = (T?W^Tc + Wq + w^q)-^ 


(19a) 

(19b) 

(19c) 


For the local model (Eq. (5)), the solution can be obtained as 
follows: 


In this derivation the response z is assumed to be deterministic. 
When the parameter uncertainties are taken into account, the cautious 
controller can be obtained by using the expected value of the per- 
formance function: 

- ® |e 

where it is assumed that W 2 is diagonal, and 0ji is deterministic. 
For the case of the open-loop control (zq feedback), there follows 

E [w2j(zjo + 

j 

L ''zj<2jo + 
j 



1 




40n»46 


A0 


n 


( 21 ) 


where 


+ D »zl(«z4 + + ejM„6„) (22) 

j 

Ptt ^tzl 

M = (23) 

L^zt ^zzj 

= (n X m) matrix 
Mtz ^zt ^ vector 

M 22 = scalar 
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So the performance function becomes 


/ \/®n\ 

J = + e>een + L + A0nWA6 (24) 

\j /M/ 


The solution for the control that minimizes J is then 


®n “ ^^0 + ^AeVi + Cq 


(25) 


The gain matrices C and C/^e are the same as for the deterministic 
controller using the identified values of the parameters and with W 0 
replaced by 


We + E W2,j(Mtt)jn 

.J 


(26) 


The new constant term is 


Co = -D E 

.j 


E W2j(Mt^) 


»] 


(27) 


Similarly, for the case of the closed-loop control (z^^_j^ feedback); 
the performance function is 


J - + E 


(28) 


The solution is identical to that for the deterministic controller 
using the identified values of the parameters and with W^e replaced 
by 

«A6 + E "zjMjn (29) 

2 

Gust Model 

In the past, studies dealing with gust-suppression systems 
have used gust models such as the von Kdrmdn model, Dryden model, and 
step, sinusoidal, and sine-squares model, etc. These models have 
some advantages as well as disadvantages. The von Kdrmdn and the 
Dryden models are expressions derived from statistical techniques in 
the frequency domain. These expressions are close to the natural 
turbulent flows in the sense of the statistics; however, they are not 
able to show the individual flow pattern of turbulence in the time 
domain. Therefore, they are not suitable for the time-wise numerical 
calculations. On the other hand, the step and sinusoidal gust 
models are very simple yet different from true gust shapes. These 
expressions are easy to handle in numerical calculations. In this 
study, the following gust representation is assumed: 

{g} = (ug* ’^G* Wq)^ = (uGo. vgo. wgo)^ exp(lo)Gt). Here uqo. vqo . 
and wqo are horizontal, lateral, and vertical gust amplitudes, 
respectively. For the step gust, oiQ - 0. For this study, the rotor 
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is assumed to have the blades at 0®, 90®, 180®, and 270® when the 
sinusoidal gust is initially encountered by the rotor (the gust veloc- 
ity field is convected past the rotor by the helicopter forward speed) . 

Numerical Calculations 

Numerical calculations have been performed to investigate this 
adaptive blade control algorithm. Blade properties and the rotor 
operating conditions are shown in Tables 1 and 2, respectively. The 
gust shape is schematically depicted in Fig. 1. The hub forces 
(Fjj» Fy* F 2 ) and the hub moments (Mjj, My, M 2 ) in the nonrotating 
frame are shown in Fig. 2. The initial estimate of the transfer 
function Tg is obtained by using a harmonic method program 
(Ref. [18]). Tables 3 and 4 show the transfer function of the rotor 
for fg = 2.0 Hz and fg = 0.0 Hz (step gust), respectively. For 
the step gust, gust frequency is assumed to be zero and the control 
inputs A0 and the outputs z are assumed as follows: 


A0 = + 0’ cos + 0g sin 

{z} = (AF 20 , AFjjo, AMyg, AFyp, AMxo» AM 20 ) 


(30) 


where the subscript 0 denotes the mean value. The transfer function 
Tq obtained by the harmonic method is regarded as the Initial esti- 
mation in the adaptive blade pitch control. In the controller, Tg is 
updated by the Kalman filter at every time-cycle. In this calcula- 
tion, it is assumed that there is no noise in the A0 measurement. 
Measurements z were contaminated by random measurement noise. Sev- 
eral sets of noise-to-signal ratios were studied before these calcu- 
lations. A noise-to-signal ratio of 0.05 was used. The initial 
values of Pg and Qg used in Eq. (19) were Pg = (9.8 x lO"*) N (or 
N-m) , Qg = (9.8 X 10^) N (or N-m) . The variance of the measurement 
noise, r, is assumed to be 9.8 N (or N-m). To keep the noise-to- 
signal ratio approximately constant, the following relationship is 
taken into account (see Ref. [ 8 ]): ~ %(3n+i/Jn)» 

’^min ’^max’ where Jn is the quadratic performance function 

at the nth time-cycle. The weighting matrices of the quadratic 
performance function are = [ 1.0 x 10 ~**], W 9 = [ 0 . 0 ], 

W^e = [1.0 X 10**]. 

Figure 3 is a schematic of the regulators used in this study: 

(a) adaptive open-loop and (b) adaptive closed-loop. In this gust 
alleviation system, the gust frequency is specified. Therefore, each 
time-cycle that updates the parameters must vary depending on the 
gust frequency. In most calculations, a 2.0 Hz gust frequency is 
used. This frequency is more than half of the rotor rotational 
speed, and somewhat larger than the typical mean atmospheric turbu- 
lence (usually below 1.0 Hz). An updated parameter estimate is per- 
formed every four revolutions of the rotor. During that time, the 
measurements z are discretely sampled. Since the measurements z 
are data in the time domain, they must be converted from time domain 
to frequency domain using the fast Fourier transform (FFT) (Ref. [19]). 
In the FFT, the resolution of the data in the frequency domain depends 
on the sampling time At (Shannon's theorem); in this calculation 
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At = 10°. Therefore, when the gust frequency £q falls between 
two points (that is, n Af < fg < (n + l)Af, n = integer, Af = frequency 
step), it is impossible to calculate the correct values by the FFT. 

In this analysis, correct values are approximated knowing the speci- 
fied gust frequency fg and using a cubic spline interpolation method. 
Furthermore, it is necessary to consider the calculation time of the 
optimal control, which determines the optimal control for the next 
time-cycle using the output data from the FFT. One time-cycle is 
8it rad (see Fig. 4). The sampling interval is (128/144) (8ir) rad. 

Data are sampled at every 10° (144 samples in four revolutions); 128 
is used in FFT sense power of 2 more efficient. Consequently the 
calculation time is assumed to be (16/144) (8tt) rad. As a result, 
there is a phase shift in the measured response at the beginning of 
the next time-cycle relative to the phase at the end of the previous 
sampling interval. This phase shift is accounted for in implementing 
the control. A global helicopter model is used. For the open-loop 
controller, both Zq and Tg are updated by the Kalman filter. Only 
Tg is updated for the closed-loop controller. 

Figures 5a to 5h show the hub gust responses with active blade 
pitch control for a vertical gust (fg = 2.0 Hz). The control law is 
deterministic, and the global helicopter model was used. Two inputs 
(0j^ and two outputs (cosine and sine elements of the thrust 

response) are considered. The thrust response (Fig. 5b) gradually 
decreases, but the modulated yawing moment response (Fig. 5g) is not 
minimized. The other responses (Fjj» Fy> My) do not show any 
influence of the controller. Here the aim of the controller is to 
reduce only the thrust fluctuation caused by the gust; it does not 
reduce the other responses. 

In Figures 6a to 6h, the hub responses are shown with active 
blade pitch control, again for the vertical gust. For these results 
the six control inputs (0^ to 0g) are used to alleviate all 12 oscil- 
latory hub responses. Measurements z involve the hub forces (F^, 

Fy, F 2 ) and the hub moments (Mx, My, Mg). Compared with Fig. 5, the 
reduction of the thrust response is significant (from 50% to 80%). 

It is observed that the fluctuation of the thrust does not reduce 
uniformly. This phenomenon can be explained as follows. The control 
inputs generated by the controller depend on the measurements. When 
the hub responses are decreased by means of the controller, the 
optimal control inputs necessarily become smaller. At some point, 
apparently, the optimal control inputs become ineffective in generat- 
ing rotor hub response. At the same time, the uncertainty of the 
identified parameters may increase because of measurement noise. In 
the yawing moment response (Fig. 6g) the convergence characteristics 
are improved relative to the two control input cases (see Fig. 5g) . 

The amplitude of the fluctuation is somewhat less than the uncontrolled 
case. Similar to the thrust response, the response sometimes diverges 
and converges. This is due to -the changes of the induced drag 
directly related to the thrust. Contrary to expectations, the 
responses (F^* Fy, Mx» My) show a slight continuous Increase in magni- 
tude. In the pitch angle change of the reference blade (Fig. 6h) , the 
change in the blade pitch with the gust frequency can be observed. 
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Figure 7 shows the time histories of the mean square thrust 
response. The solid line denotes the response without control. The 
dashed line corresponds to Fig. 5 (2 inputs and 2 outputs) and the 
broken line corresponds to Fig. 6 (6 inputs and 12 outputs). -As 
explained before, the parameters are updated at every four revolu- 
tions. During one time-cycle, the control inputs are kept constant. 

The response for the case of thrust control only is at first con- 
stant, and then gradually decreases. The response is reduced by 
almost 50% after 30 revolutions (five time-cycles). The response for 
the case of complete hub response control decreases to a much lower 
value, although not continuously. A 50% to 80% reduction of the thrust 
response is attained for the case of complete hub response control. 

From these results, it is found that the measurements z should Include 
not only the thrust response but also the other hub responses. 

Figures 8 a to 8 h show the hub gust responses for a vertical 
step gus£. In this case, the adaptive closed-loop controller is used. 
The measurements z involve only the longitudinal response of the 
rotor; that is, the horizontal force (Fj^) , the thrust (F 2 ) and the 
pitching moment (My). The control inputs have three elements (6{, 

83 ) according to Eq. (30). All parameters (zq and Tq) are updated at 
every rotor revolution by the Kalman filter. Significant reduction 
(nearly 100%) in the thrust response (Fig. 8 b) is achieved by the 
controller. However, compared with the uncontrolled case, the hori- 
zontal and side force responses and the rolling and pitching moment 
responses transfer to slightly different steady state conditions. In 
the yawing moment response (Fig. 8 g) , there is a significant change 
in the steady state condition. This phenomenon comes from the strong 
effect of the controller on the gust alleviation system. The con- 
troller generates the optimal control inputs in order to reduce the 
fluctuation from the steady value. If these control inputs are sig- 
nificant values relative to the baseline trim pitch inputs, a new 
trim condition could result. Referring to Fig. 8 h, the maximum 
change of the pitch angle is about 2®. The cyclic pitch angle for 
the trim condition is less than 2®. Therefore, the pitch angle change 
by the controller has an effect on the rotor trim condition. The 
flight trim condition changes and transfers to the new steady state 
condition. 

In Fig. 9, the sensitivity of the adaptive closed-loop con- 
troller performance to the transfer function is investigated for the 
step gust. Arbitrary small Initial values for the transfer function 
are used in this calculation. Again, only longitudinal hub responses 
are included in the plant model. Compared with Fig. 8 , the thrust 
fluctuation (Fig. 9b) decreases very slowly. The other responses, 
except yawing moment, remain at the same trim condition. In the case 
of the yawing moment (Fig. 9g) , the convergence to a new steady state 
can be seen after showing a significant transient. This phenomenon 
is due to the Induced drag, as mentioned before. From these results, 
it can be concluded that the performance of the controller depends on 
the initial estimate of the rotor transfer function. If a more accu- 
rate transfer function is used, quick convergence of the responses 
can be obtained. The transfer function derived by the harmonic method 
provides a reasonable initial estimate. 
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Figures 10a to lOh show the hub responses for the vertical step 
gust with the adaptive open-loop controller. For these results, the 
measurements z involve all six responses of the rotor. All elements 
of Tq are updated at every rotor revolution. The magnitudes of the 
controller Inputs are constrained by a prescribed maximum value, 

A0max “ 2.0°. For the thrust response (Fig. 10b), the deviation from 
the steady value of the thrust only gradually decreases after respond- 
ing to the step gust. The other hub forces together with the rolling 
and pitching moment remain steady. The yawing response deviates 
slowly from the steady value because of the thrust response. Compared 
with the case of Fig. 8, the pitch angle change is very slow and 
small. For the open-loop controller, the optimal control inputs 
depend on the uncontrolled response Zq. If the Zq are small, then 
the optimal control Inputs may be small. For this case, the optimal 
control Inputs at the beginning stage of controller operation are 
less than 0.5°. 

In Fig. 11, the time histories of the quadratic performance 
function J are shown for the vertical step gust. A solid line 
denotes the time history of J for the uncontrolled case. The J 
with control are shown for both the adaptive closed-loop control 
(corresponding to Fig. 8) and the adaptive open-loop control (corre- 
sponding to Fig. 10). For the adaptive closed-loop controller, the 
J decreases rapidly after showing a sharp Increase. The J for the 
adaptive open-loop controller decreases moderately after showing the 
same increase. Both control schemes demonstrate the effect the con- 
troller has on the reduction of the hub gust responses. 

For the deterministic controller considered in this analysis, 
the uncertainty of the parameter identification has not been accounted 
for in the calculation of the optimal control inputs. The cautious 
properties can play an Important role in helicopter vibration reduc- 
tion schemes (Ref. [14]). In Figs. 12a to 12h, the hub gust responses 
with the cautious controller are shown for the vertical sinusoidal 
gust (fg = 2.0 Hz). Referring to the results with the deterministic 
controller (Fig. 6), there is no significant difference between these 
controllers. The cautious controller involves the variance of the 
error in the estimate of the transfer function T before the measure- 
ment. The more uncertain the parameter estimates, the larger the 
variance of the error. In this case, the term concerned with the 
cautious properties 

J 

in Eq. (26) was about 1% of the term T^Wglg in Eq. (19c). Therefore, 
the cautious properties have little effect on the controller perfor- 
mance in this study. 

Results from applying adaptive blade pitch control for a three- 
dimensional gust are shown in Figs. 13a to 13h. The three gust com- 
ponents have amplitudes Ugg =3.0 m/sec, vqq = 2.0 m/sec, and 
Wqq = 1.8 m/sec; all have frequency fg = 2.0 Hz. A large reduction 
in the thrust response is observed. In general, all the hub responses 
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are similar to the case of a vertical gust only (Fig. 6). In this 
study it was observed that the rotor hub responses were most sensitive 
to the vertical gust component. For the three component gust, the 
controller's cautious properties are not considered because of their 
demonstrated minimal effect on controller performance for the one 
component gust case. 


3. Conclusions 


Feedback control systems to reduce the gust-induced vibration 
of a helicopter rotor have been studied by an adaptive blade pitch 
control. For the three gust components, the horizontal and lateral 
gust components have little influence on the hub gust response. How- 
ever, the vertical component of gust has a significant Influence. In 
studying the gust response the rotor thrust shows the most significant 
change compared with other hub responses. Hence, reducing the thrust 
response may be the only necessary aim of the gust alleviation system. 
In this study, the major, the subharmonic, and the superharmonic 
inputs in the nonrotating frame are considered in the controller. As 
the gust frequency Increases, higher frequency terms of the gust 
response become large. To alleviate these high frequency terms in 
the vibratory response, the number of terms included in the control 
Inputs increases, making the controller more complex. The gust fre- 
quency is prescribed at the Initial stage of the calculations. There- 
fore, this type of regulator does not apply to random gust responses. 
From these theoretical results, the following conclusions can be drawn. 

1. The performance of the controller depends on the initial 
estimate of the rotor transfer function. When the exact transfer 
function is used, the convergence of the responses can be achieved 
quickly. The transfer function derived by the harmonic method 
(Ref. [18]) gives a good estimate for initializing the controller. 

2. In using the FFT to convert measurements from the time 
domain to the frequency domain, some approximation must be made because 
of the gust frequency. This may increase the uncertainty of the mea- 
surement. A more accurate frequency domain determination method is 
required. 


3. Using a cautious controller to account for the uncertainty 
of the parameter identification has little influence on the improve- 
ment of the regulator performance. 

4. For the case of a sinusoidal gust, the adaptive open-loop 

regulator is best suited for the gust alleviation system. Results 
show that a 50% to 80% reduction of the thrust response can be obtained. 
The regulator studied in this paper is shown to be applicable to a 
three dimensional gust. ' 

5. For the step gust, the adaptive closed-loop regulator per- 
forms better than the adaptive open-loop regulator. The closed-loop 
regulator yields a rapid reduction of the gust-induced thrust response 
(almost 100%), even though it violates the trim condition. The open- 
loop regulator shows that convergence of the thrust response is slow. 
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TABLE 1.- ROTOR CHARACTERISTICS 


Rotor radius, R 
Blade lift slope, a 
Blade semichord, b 
Number of blades, N 
Rotor rotational speed, 
Blade twist rate, 0^. 
Position of flapping hinge. 
Blade cutoff, r^^ 

Blade c.g. position, r^g 
Position of lag damper, r^ 
Lag damper coefficient, C^ 
Blade mass, m 

Moment of inertia of blade. 
Mass moment of blade. Mg 
Lock number, y 
Wing section 
Gross weight, W 


8.53 m 

5.73 

0.2185 m 
4 

23.67 rad/sec 
- 8 ° 

rg 0.3 m 
0.594 m 

2.74 m 
0.3 m 

1000.0 N-m-sec/rad 

106.4 kg 

Ig 1593.9 kg-m^ 

1659.1 kg-m^ 

8.84 

NACA 0012 

62259.4 N 
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TABLE 3 


[Tc] = 


TABLE 2.- OPERATING CONDITIONS 


Advance ratio, ]i 0.18 

Collective pitch angle, 0 q 8.0° 

Longitudinal pitch angle, 6j^g -1.48° 
Lateral pitch angle, 0.695° 

Inflow ratio, X 0.0179 


.- TRANSFER FUNCTION OF ROTOR FOR GUST FREQUENCY OF 2. 


-25.52 

-576.5 

576.5 

-25.52 

9.831 

-59.62 

66.51 

15.42 

34.46 

-112.7 

112.7 

34.46 

-63.17 

2.580 

-2.580 

-63.17 

-1185 

4.166 

-4.166 

-1185 

969.5 

77.35 

.-77.35 

969.5 


-192.6 

22.40 

22.40 

-192.6 

-211.6 

-44.92 

3.635 

-179.3 

2397 

-85.89 

85.89 

2397 

-48.63 

216.3 

-216.3 

-48.63 

96.97 

2508 

-2508 

96.97 

-10.06 

-293.7 

293.7 

-10.06 


81.50 

-11.92 

11.92 

81.50 

50.86 

53.35 

-88.02 

59.16 

2561 

91.49 

-91.49 

2561 

-56.87 

37.14 

-37.15 

-56.87 

87.10 

-2618 

2618 

87.10 

4.302 

202.5 

-202.6 

4.305 


TABLE 4.- TRANSFER FUNCTION OF ROTOR 
FOR STEP GUST 


"0.996 

8.545 

2.984 

-0.496 

8.004 

2.077 

1.849 

1.872 

0.635 

68.42 

252.3 

95.31 

-16.62 

289.0 

81.92 

.8.599 

17.79 

5.652 


0 Hz 


X 10^ 
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COORDINATE SYSTEMS 


ROTOR: Xr, Yr, Zr 
HUB:Xh.Yh.Zh 


^R'^H 



Fig. 1, Gust shape and coordinate systems. 





Fig. 2. Hub forces and moments. 
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(a) Open-loop regulator. 
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(b) Closed-loop regulator. 

Fig. 3. Adaptive regulators for gust alleviation 
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Fig. 4. Sequence of controller Implementation. 
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Fig. 6. Time history of hub responses with adaptive open-loop blade 
pitch control for a vertical sinusoidal gust: wqq = 1.8 ra/sec, 

fQ = 2.0 Hz. 
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Fig. 7. Time history of the mean square thrust. 
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Fig. 9. Time history of hub responses with adaptive closed-loop blade 
pitch control for a vertical step gust; wnn =1.8 m/sec. 
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Fig. 10. Time history of hub responses with adaptive open-loop blade 
pitch control for a vertical step gust: wqq =1.8 m/sec. 
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WITHOUT CONTROL 



Fig. 11. Time history of quadratic performance function J. 
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Fig. 12. Time history of hub responses with adaptive open-loop blade 
pitch control with caution for a vertical sinusoidal gust: 

Wqp = 1.8 m/sec, Iq = 2.0 Hz. 
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